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It had been conjectured that Diffusion Limited Aggregates and Laplacian Growth Patterns (with 
small surface tension) are in the same universality class. Using iterated conformal maps we construct 
a 1-parameter family of fractal growth patterns with a continuously varying fractal dimension. This 
family can be used to bound the dimension of Laplacian Growth Patterns from below. The bound 
value is higher than the dimension of Diffusion Limited Aggregates, showing that the two problems 
belong to two different universality classes. 



Laplacian Growth Patterns are obtained when the 
boundary F of a 2-diniensional domain is grown at a 
rate proportional to the gradient of a Laplacian field P. 
Outside the domain V^P = 0, and each point of F is ad- 
vanced at a rate proportional to VP It is well known 
that without ultra-violet regularization such growth re- 
sults in finite time singularities In correspondence 
with experiments on viscous fingering one usually adds 
surface tension, or in other words solves the above prob- 
lem with the boundary condition P = an where a is the 
surface tension and k the local curvature of F The 
other boundary condition is that as r — > oo the flux is 
VP — const xf/r. Fig. 1 (left) shows a typical Laplacian 
Growth Pattern. 




FIG. 1. Left: Typical Laplacian Growth pattern with sur- 
face tension, cf. Ref.[6]. Right: Typical DLA cluster of 
100000 particles. 



Diffusion Limited Aggregation (DLA) |^ begins with 
fixing one particle at the center of coordinates in 2- 
dimensions, and follows the creation of a cluster by re- 
leasing fixed size random walkers from infinity, allowing 
them to walk around until they hit any particle belong- 
ing to the cluster. Since the particles are released one by 
one and may take arbitrarily long time to hit the cluster, 
the probability field is quasi-stationary and in the com- 
plement of the cluster we have again V^P — 0. In this 
case the boundary condition on the cluster is P = 0, but 
finite time singularities are avoided by having finite size 



particles. The boundary condition at infinity is exactly 
as above. A typical DLA is shown on Fig.l (right). 

In spite of the different ultra-violet regularizations of 
Laplacian Growth and DLA, it was speculated by many 
authors |^ that the two problems belong to the same 
"universality class" , and it was expected that the result- 
ing fractal patterns will have the same dimension. In this 
Letter we argue that this is not the case: there are deep 
differences between the two problems, and in particular 
Laplacian Growth Patterns have a dimension consider- 
ably higher than DLA. In one sentence, the differences 
between the problems stem from the fact that Lapla- 
cian Patterns are grown layer by layer, whereas DLA 
is grown particle by particle. Unfortunately, traditional 
techniques used to grow Laplacian Growth patterns, ei- 
ther numerical |^ or experimental Q , fail to achieve pat- 
terns large enough to extract reliable dimensions (and 
see Fig. 9 in ^ for example). The numerical algorithms 
are extremely time consuming due to the stiffness of 
the equations involved; experimentally it is difficult to 
construct large quasi two-dimensional (Hele-Shaw) cells 
without introducing serious deformations. 

The aim of this Letter is to provide a scheme to simu- 
late the zero surface tension Laplacian Growth that has a 
finite size regularization and thus does not suffer from fi- 
nite time singularities. We introduce a 1-parameter fam- 
ily of growth processes based on iterated conformal maps 
Contrary to DLA which grows particle by parti- 
cle, we will construct the family of growth processes to 
mimic Laplacian Growth, in which a layer is added to the 
boundary F at each growth step, with a width propor- 
tional to the gradient of the field. Consider then $(")(k;) 
which conformally maps the exterior of the unit circle e'^ 
in the mathematical w-plane onto the complement of the 
(simply-connected) cluster of n particles in the physical 
^-plane. The unit circle is mapped onto the boundary of 
the cluster. The map $("^(w) is made from compositions 
of elementary maps (px.e, 

$(")(«;) =$(«-l)(0;,^ ,^(u;)) , (1) 

where the elementary map 4>\,0 transforms the unit circle 
to a circle with a "bump" of linear size -s/A around the 
point w = e'^. In this Letter we employ the elementary 
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map 



-(l + w) 




With this choice the map 

$(")(^^^) adds on a new semi- 
circular bump to the image of the unit circle under 
$("-!) (y;), xhe bumps in the x-plane simulate the ac- 
creted particles in the physical space formulation of the 
growth process. The recursive dynamics can be repre- 
sented as iterations of the map (j)\^_g^{w)^ 



(4) 



With the present technique it is also straightforward 
to determine the dimension. The conformal map 

$(»)(^) 

admits a Laurent expansion 



(n) 
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The coefficient of the linear term is the Laplace radius, 
and was shown to scale like 



(") 
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where S is the area of the cluster (the sum of the actual 
areas of the bumps in the physical space). On the other 

hand f[^^ is given analytically by 



Fr = n votm , 



(7) 



fc=i 



and therefore can be determined very accurately. 

Different growth processes can be constructed by 
proper choices of the itineraries jl^l, and rules 

for determining the areas of the bumps {A^J^^j^. In DLA 
growth 1^,^ one wants to have fixed size bumps in the 
physical space, say of fixed area Aq. Then one chooses in 
the nth step 
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Ao 



|$(n-l)'(gie„)|i 



DLA growth 



The probability to add a particle to the boundary of 
the DLA cluster is the harmonic measure, which is uni- 
form on the circle. Thus in DLA the itinerary {Oi}^^^ is 
random, with uniform probability for 9i in the interval 
[0,2^]. 

For our present purposes we want to grow a layer of 
particles of varying sizes, proportional to the gradient of 
the field, rather than one particle of fixed size. This en- 
tails three major changes. First, if we want to grow one 
particle of size proportional to the gradient of the field, 



(i.e. area proportional to j^*^" '^)'(e*^")| ^) we need to 
choose 



An — 



Ao 



|$(n-l)'(g»e,.)|4 



present models . (9) 



Second, to grow a layer, we need to accrete many particles 
without updating the conformal map. In other words, to 
add a new layer of p particles when the cluster contains 
m particles, we need to choose p angles on the unit circle 
{(^m+k}k=i- At these angles we grow bumps which in the 
physical space are proportional in size to the gradient of 
the field around the m-particle cluster: 



|$(m)'(gie„+fe)|4 



1,2. 



(10) 



Lastly, and very importantly, we need to choose the 
itinerary {^m+fe}fc=i which defines the layer. This 
itinerary is chosen to achieve a uniform coverage of the 
unit circle before any growth takes place. The parame- 
ter that will distinguish one growth model from another, 
giving us a 1-parameter control, is the degree of coverage. 
In other words we introduce the parameter 
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fc+i 



This parameter is the fraction of the unit circle which is 
covered in each layer, with the limit of Laplacian Growth 
obtained with C = 1. It turns out to be rather time con- 
suming to grow fractal patterns with C close to unity. 
But we will show below that this is hardly necessary; al- 
ready for C of the order of 1 /2 we will find patterns whose 
fractal dimension significantly exceeds that of DLA, of- 
fering a clear lower bound on the dimension of Laplacian 
Growth patterns. 

Once a layer with coverage C had been grown, the field 
is updated. To do this, we define a series {Ok}^^i accord- 
ing to 



(12) 



Next we define the conformal map used in the next layer 
growth according to 



(8) $(™+P)(tj) = $(™)o0e„ 



■ ° 9e„ 



(13) 



It is important to notice that on the face of it this con- 
formal map appears very similar to the one obtained in 
DLA, Eqs.(0), But this is deceptive; the distribution 
of 9 values is different, we do not update the map after 
each particle, and the growth rule is different. 

We can achieve a uniform coverage C using various 
itineraries. One way is to construct the "golden mean 
trajectory" 9m+k+i = Om+k + '^T^P where p = (V5- l)/2. 
At each step we check whether the newly grown bump 
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may overlap a previous one in the layer. If it does, this 
growth step is skipped and the orbit continues until a 
fraction C is covered. Another method is random choices 
of 0,„+fe with the same rule of skipping overlaps. We 
have tried several other itineraries. Of course, to be an 
acceptable model of Laplacian growth the resulting clus- 
ter should be invariant to the itinerary. This invariance 
is demonstrated below. The central thesis of this work 
is that the dimension of the resulting growth patterns 
is dependent on C only, and not on the itinerary chosen 
to achieve it. Numerically it is more efhcicnt to use the 
golden mean itinerary since it avoids as much as possible 
previously visited regions. In order to achieve compa- 
rablegrowth rates for different layers we inflated Ao in 
Eq. ( |lO[) according to Aq — )■ mAo in the layer composed of 
p particles {m + k}^.^^. In Fig. 2 we show Fi of clusters 
grown by choosing 3 different itineraries to produce the 
layers and for two values of C. 
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FIG. 2. Log-log plots of Fi vs. S of six individual clus- 
ters, using 3 differnt itineraries for layer construction, with 
two values of C C = 0.3 (upper group) and C — 0.5 (lower 
group). Here we use the golden-mean, random and the period 
doubling itineraries (see Ref.[10]). 



We conclude that the dimension (determined by the 
asymptotic behavior of Fi vs. S) does not depend on the 
itinerary used to form the layers but on C only. 



In Fig. 3 we show three fractal patterns grown with this 
method, with three different values of C. Even a cursory 
observation should convince the reader that the dimen- 
sion of these patterns grows upon increasing C. 





FIG. 3. Patterns grown with 3 different values of C by us- 
ing the golden-mean itinerary: a) C = 0.1, b) C = 0.3, c) 
C = 0.5. 



In order to calculate the dimension we averaged F\ of 
many clusters produced by the golden mean itinerary, 
each with another random initial angle in each layer. 
Plots of the averages (F\) for 3 values of C are presented 
in Fig. 4. 
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FIG. 4. Linear regressions of log-log plots of (Fi) vs. S for 
3 values of C: 0.1 (solid line), 0.3 (dotted) and 0.6 (dashed). 
The slopes of the curves imply dimensions D=1.37, D=1.75 
and D=1.85 respectively. The averages are taken over at least 
20 clusters. 
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We conclude that the dimension of the growth pat- 
tern increases monotonically with C, with D « 1.85 when 
C = 0.6. 

The main point of this analysis is that the dimension 
of Laplacian growth patterns is bounded from below by 
the supremum on the dimensions obtained in this fam- 
ily of models. First, Laplacian Growth calls for C — 1. 
Second, in Laplacian Growth the boundary condition is 
P ~ an, suppressing growth at the tips (and relatively 
favoring growth in the fjords) compared to growth with 
the boundary condition P = 0. Accordingly, on the basis 
of the results shown in Fig. 4, we propose that the dimen- 
sion of Laplacian Growth patterns exceeds 1.85, putting 
it distinctively away from the dimension of DLA which 
is about 1.71 Q. 

In hindsight, it is difficult to understand how the con- 
sensus formed in favor of DLA and Laplacian Growth 
being in the same universality class. Superficially one 
could say that in DLA the update of the harmonic mea- 
sure after each particle is not so crucial, since the effect 
of such an update is relatively local [ p2[ . Thus it may 
just work that a full layer of particles would be added 
to the cluster before major interaction between different 
growth events takes place. However this view is com- 
pletely wrong. An incoming random walker lands on top 
of a previously attached one very often. To see this, con- 
sider how many angels {Oj} can be chosen randomly on 
the unit circle before the first overlap between bumps (of 
linear sizes = ■y/A„(e*^j)). To get the order of magni- 
tude take ej — e — {^/Xn). The average number of times 
that we can choose randomly an angle before the first 
overlap is M{e) ^ The Length of the unit circle that 
is covered at that time by the already chosen bumps is 
C{e) = eA/'(e) It was shown in B that for DLA 

(A„) ^ i, so that e ~ implying N(n) ~ n^/^. No- 
tice that this result means in particular that for a DLA 
cluster of 1 million particles only less than 50 random 
walkers can be attached before two of them will arrive at 
the same site! Moreover, C{n) ^ ^ for n ^ oo, 
which means that as the DLA cluster grows, our cover- 
age parameter C goes to zero, rather than to unity where 
Laplacian Growth is. Taking spatial fluctuations of A„ 
into account may change the exact exponents but not the 
qualitative result. This argument clarifies the profound 
difference between growing a whole layer simultaneously 
and particle-by-particle. Note however that DLA is NOT 
the C — > limit of our 1-parameter family due to the dif- 
ference between Eqs.(||) and (^. 

The results of this study underline once more the deli- 
cacy of the issues involved. Fractal patterns depend sen- 
sitively on the details of the growth rules. Even though 
the analytic presentation seems very similar, to the de- 
gree that many researchers were led to believe in wide 



universality classes, we showed here that one must be 
much more cautious. By lifting the models into fam- 
ilies of growth patterns depending on a parameter we 
could demonstrate strong variability of the fractal dimen- 
sion. Here we constructed the family to bound from be- 
low Laplacian Growth patterns. A similar family can be 
constructed to bound DLA from above. This and other 
aspects of this method will be reported elsewhere. 
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